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1.  INTRODUCTION 


The  objective  of  this  study  progrcun  has  been  to 
study  methods  of  detection  of  coherent  or  partially 
coherent  optical  radiation  in  a noncoherent  background. 
This  study  is  primarily  intended  for  application  at 
low  light  levels  and  low  ratios  of  coherent  to  non- 
coherent (i.e./  signal/noise)  light.  This  study  has 
application  for  identifying/locating  a source  of 
partially  coherent  light  (e.g.,  laser)  in  daylight  or 
in  the  presence  of  thermal  background  radiation. 

Essentially  this  measurement  method  is  based  upon 
measurements  of  temporal  coherence  of  received  light. 
The  background  light  is  almost  totally  noncoherent. 
Laser  radiation  is  partially  coherent,  and  can,  in 
certain  cases  approach  total  coherence. 
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2.  REVIEW  OF  METHODS  FOR  MEASURING 


TEMPORAL  COHERENCE 

One  of  the  classical  methods  of  measuring  temporal 
coherence  is  based  upon  a form  of  interferometry.  In 
such  methods  the  beam  of  light  whose  coherence  is  to 
be  measured  is  separated  into  two  beams.  One  of  the 
beams  is  time  delayed  relative  to  the  other  (e.g.,  by 
introducing  a path  length  difference) , and  then  the  becuns 
are  recombined  in  such  a way  as  to  form  visible  interference 
fringes.  The  interference  fringe  contrast  ratio  is  a 
function  of  the  temporal  coherence  which  is  being  sought. 

This  interferometric  measurement  could,  in  principle, 
be  used  as  a method  of  detecting  laser  radiation.  However, 
it  tends  to  be  useful  for  detecting  a temporally  coherent 
source  only  at  large  signal  levels  and  large  signal/noise 
compared  to  the  intent  of  our  investigation.  Thus  our 
study  will  be  devoted  more  to  techniques  of  measuring 
temporal  coherence  which  operate  at  the  lowest  possible 
levels. 

There  are  several  methods  for  measuring  temporal 
coherence  which  are  applicable  to  low  signal  levels. 

Each  of  these  involves  measurements  of  statistical 
parameters  of  the  photon  arrival  rate.  The  signal 
detection  method  which  we  have  studied  in  our  investigation 
uses  photon  arrival  rate  statistics  to  estimate  various 
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moments  of  the  statistical  distribution  function  for  the 
photon  arrival  rate.  We  will  demonstrate  shortly  that 
these  are  distinctly  different  for  coherent  and  noncoherent 
radiation . 

In  a later  section  of  this  report  we  describe 
instrumentation  and  procedures  for  estimating  the  moments 
of  the  photon  arrival  rate  distribution  function  from 
finite  numbers  of  observations.  From  these  we  describe 
a signal  processing  method  which  leads  to  a definition 
of  a signal-to-noise  ratio  from  which  the  signal  detecta- 
bility can  be  assessed.  We  further  determine  the  size  of 
the  Scimple  space  required  for  acceptable  confidence  in 
the  detection  of  the  coherent  radiation. 


3.  BRIEF  TUTORIAL  REVIEW  OF 


TEMPORAL  COHERENCE 

The  temporal  coherence  of  an  optical  field  is 
a measure  of  its  ability  to  produce  stable  time 
independent  interference  fringes.  The  temporal  coherence 
can  equally  well  be  represented  in  the  time  or  frequency 
domain  using  classical  statistical  concepts. 

In  classical  statistics  the  temporal  coherence  of 
a scalar  stationary  random  process  V(t)  is  given  by 
the  autocorrelation  function  Ry('^)  or  the  power  spectral 
density  W^(w) 

R„(T)  « Lim  i / V(t)  V(t  + T)dt 

''  T-»-«>  ^ *'0 

Wy(u)  = y*  Ry(T)  e~^“'^  dT 

In  the  semi-classical  approximation  (i.e.,  ignoring 
quantum  effects)  these  same  functions  are  useful  for 
representing  the  temporal  coherence  of  any  Cartesian 
component  of  an  optical  electric  or  magnetic  field. 

A linearly  polarized  optical  field  is  conveniently 
represented  by  the  complex  amplitude  U which  is 
related  to  the  intensity  I by; 

I = |ur  = uu 
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This  scalar  amplitude  is  proportional  to  the  electric 
field  peak  amplitude  E 


U = 


"0 


where 


E = ||E|| 

E = vector  electric  field 


Hg  = wave  impedance 


In  the  above  both  U and  E are  complex  analytic  functions 
which  are  derived  from  the  corresponding  real  electric  fields 
The  coherence  of  an  optical  field  is  represented 
by  the  mutual  coherence  function  T(r^,  ?2»  t)  where: 

T 

r(r^,  r^,  T)  = Lim  ^ J" U*(rj^,t)  U(?2»  t + t)dt 

T-^oc  - _ 


where 

r.  « = vector  positions 

X / ^ 

T * time  difference 

It  is  convenient  to  denote  r(r^,  r2^  x)  by  the  symbol 
this  entire  discussion  we  are  assuming 
that  U is  ergodic. 

Prom  this  definition  of  the  mutual  coherence  function 
the  corresponding  frequency  domain  function  GCr^^,  r2,  v) 
can  be  found: 


This  function  is  called  the  mutual  spectral  density  or 
sometimes  the  cross  spectral  density. 

In  most  calculations  it  is  helpful  to  use  normalized 
functions.  The  normalized  mutual  coherence  function  is 
called  the  degree  of  coherence  and  is  denoted 
Y (r^r2T)  = Yj^2 


r (r-r^T) 

{rCrj^r^^O)  rr2r20}-^/^ 

Note  that  bounded; 

The  corresponding  normalized  spectral  density  g^^2 
is  given  by 
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(V) 


= /: 


Yjl2 


-i2nvT 


dT 


(V) 


There  is  an  important  special  case  when  considering 
light  from  a laser.  Laser  generated  light  occupies  a 
very  narrow  spectral  width  Av.  In  this  case  the 
normalized  spectral  density  can  be  approximated  by: 
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E = entire  source  area  which  is 
radiating  light 

This  particularly  important  special  case  is  known 
as  the  Van-Citterr  - Zernike  theorem.  It  has  the  importance 
of  separating  the  spatial  coherence 

temporal  coherence  which  is  given  by  In  the 

remainder  of  our  work  we  will  make  the  simplifying 
assumption  that  the  light  whose  temporal  coherence  we 
are  attempting  to  measure  occupies  a very  narrow  spectrum 
Av  compared  to  the  average  frequency  <v>.  Thus  in 
this  study  we  are  attempting  to  measure  by 

various  methods  rather  than  attempting  to  find 

The  method  which  we  have  been  studying  for  the 

measurement  of  y-.i  (i)  is  based  upon  measurements  of 

xl 

the  statistics  of  photon  arrival  rates.  This  study 
has  been  motivated  by  the  work  in  Ref.  4 which  showed 
the  relationship  between  and  photon  statistical 

measurements.  There  it  is  shown  that  for  light  having 
Gaussian  amplitude  statics  we  have; 

Pc(t)  = n < I(t)  >[!■'■  1^11^'^^  1^] 

where 

P(,(t)  = probability  of  receiving  a 

photon  at  time  t + t given 
a photon  at  time  t 


JL 


utetfliHitaMk  - 


- ^ 


<I(t)>  = average  light  intensity 

n = detector  efficiency 
<n>  = average  number  of  photons 

generated/ unit  time  by  <I(t)> 

This  formulation  is  not  strictly  applicable  to  the 
detection  of  laser  light  because  the  latter  does  not 
have  Gaussian  amplitude  statistics.  However  it  is  shown 
(Ref.  4)  that  noncoherently  generated  photons  have  a 
tendency  to  bunch  together  having  a nonuniform  arrival 
rate;  whereas  perfectly  coherent  constant  amplitude  light 
has  a uniform  arrival  rate.  This  result  suggests  that 
the  arrival  rate  statistics  are  different  for  coherently 
generated  photons  that  for  noncoherently  generated  photons 
A significant  portion  of  the  research  conducted  during 
this  study  program  has  been  devoted  to  a careful  investi- 
gation of  the  statistics  of  photon  arrival  rates. 

Section  5 of  this  report  describes  the  results  of 
our  theoretical  study  of  photon  arrival  rate  statistics. 
There  it  is  shown  that  certain  of  the  statistical 
moments  for  the  photon  arrival  rate  are  sufficiently 
different  for  coherent  and  noncoherent  light  to  provide 
the  basis  of  detection  of  the  former  when  concealed  by 
the  latter.  In  Section  6 we  define  a signal  quantum 
from  which  it  is  possible  to  detect  coherent  light  in  a 
noncoherent  background. 


4.  INSTRUMENTATION  FOR  OPTICAL  SIGNAL 


DETECTION  METHOD 


The  method  which  we  are  studying  for  detection 
of  coherent  radiation  can  be  understood  with  reference 
to  Fig.  1. 


Statistical 

Parameters 


Figure  1 

* *1 

This  figure  is  a block  diagram  of  a portion  of  the 
optical  signal  detection  instrumentation.  This 
instrviment/  which  is  presumed  to  be  receiving  light  from 
some  region  of  space/  receives  light  having  intensity 
I^^^  from  that  space.  This  region  of  space  is  that 
portion  contained  inside  a solid  angle  which  is  deter- 
mined by  the  instrument  optics. 

The  optics  include  a lens  system  which  might, 
for  example,  be  a telescope  and  possibly  some  sort  of 
preselector  or  filter  such  as  a Fabry  Perot  etalon. 

The  parameters  of  these  optical  components  are  a function 
of  the  intended  application  of  the  system.  For  example. 
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be  reasonable  to  remove  the  light  outside  of  the 
inhomogeneously  broadened  linewidth  for  that  particular 
laser  and  wavelength  by  means  of  an  optical  filter. 

Referring  to  the  previous  section  we  assumed 
that  the  spectral  width  of  the  optical  source  Av  is 
small  compared  to  the  mean  frequency  <v>.  The  preselector 
filter  which  we  mentioned  above  is  a means  of  achieving 
this  narrow  spectral  profile. 

Other  factors  which  influence  the  specification  of  the 
parameters  of  the  system  optics  are;  the  desired  region 
of  coverage,  the  total  number  of  wavelengths  which  must  be 
covered  and  the  corresponding  linewidths/bandwidths.  In 
this  study  we  presume  that  the  incident  light  is  filtered 
about  a center  frequency  where  v is  the 

d> 

average  optical  frequency  over  the  filter  passband 
(i.e. , v_  = <v>) . 

a 

The  light  which  emerges  from  the  system  optics 
reaches  the  active  area  of  a photo  detector.  In  our 
study  we  are  considering  low  light  levels  and  low 
signal/noise.  A detector  of  great  sensitivity  will  be 
required  for  this  purpose  which  in  limiting  cases 
should  be  capable  of  detecting  individual  photons. 

That  is,  the  detector  output  signal  must  rise  above  the 
quiescent  noise  level  in  the  presence  of  a single  photon. 

Moreover  we  will  be  concerned  with  rapid  variations 
in  the  intensity  of  the  light.  The  detector  should  be 
capable  of  responding  to  variations  in  light  intensity 


which  occur  in  intervals  on  the  order  of  nanoseconds. 


This  implies  a detector  having  a bandwidth  of  the 
order  of  a few  GHz.  The  combination  of  large  effective 
detector  bandwidth  and  sensitivity  suggests  that  the 
detector  will  probably  have  to  be  a cooled  photo  multiplier, 
However,  it  is  not  the  purpose  of  this  investi- 
gation to  develop  a design  for  an  instrument.  Therefore 
we  will  not  pursue  the  characteristics  of  the  detector 
further  than  to  say  that  it  must  be  capable  of  detecting 
very  small  numbers  of  photons  and  must  have  a bandwidth 


in  excess  of  about  10  Hz. 

In  the  remainder  of  this  report  it  will  be  assumed 
that  the  detector  is  a photo  multiplier  having  a 
multiplication  factor  That  is  the  ratio  of  anode 

current  i^  to  cathode  current  is  given  by: 


-rii  = M 
Ic 


The  cathode  current  is  determined  by  the  number  of 
incident  photons,  by  the  quantum  efficiency  of  the 
cathode  and  by  the  optical  frequency.  Assuming  that  the 
average  energy  per  photon  in  the  filtered  light  reaching 
the  cathode  (i.e.,  hv  ) is  sufficient  to  overcome  the 

a 

cathode  work  function,  then  some  nonzero  fraction  of 
the  incident  photons  will  generate  free  electrons  at 
the  cathode.  Denoting  the  number  of  incident  photons/ 


unit  time  which  strike  the  cathode  N and  the  number 

P 

of  liberated  electrons/unit  time  N^,  then 
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N 


e 


where  C is  the  cathode  quantum  efficiency.  It  is 

convenient  to  call  N the  photon  flux  rate. 

P 

This  photon  flux  rate  is  linearly  proportional 
to  intensity  I(t): 


I ds 


where  is  the  active  area  of  the  detector  which 

received  intensity  I.  In  v/riting  this  expression  we 
assume  normal  incidence  upon  a planar  cathode.  For 
our  analysis  we  make  the  simplifying  assumption  that 
I is  uniform  over  for  which  v/e  obtain: 


The  cathode  current  is  i_  given  by 


i = g N 
c ^ e 


auid  the  anode  current  is  given  by; 


i = q N M 
a e pm 


In  writing  these  expressions  we  have  ignored  the 
transit  time  of  electrons.  Ignoring  these  effects  is 
equivalent  to  assuming  that  the  detector  has  infinite 
detection  bandwidth.  We  will  make  this  assumption 
throughout  this  report. 


14 


Our  signal  detection  method  is  based  upon  the 
statistics  of  the  photon  flux  rate.  In  collecting  data 
for  the  required  statistical  calculations  it  is  necessary 
to  sample  the  detector  output  current.  This  is 
accomplished  by  means  of  a very  fast  electronic  switch. 
Denoting  the  switched  output  current  i^ 
the  mth  sample  we  can  write 


S/m 


(t) 


i,  (t)  t < t < t + T 

cl  m — m 


= 0 


t + T < t < 
m — 


^m+1 


t * mr 
m 

m ^ 1/  2/  ..#/  M 

T = sampling  period 
T = sampling  interval 
T > T 

M = number  of  samples  taken 


It  will  be  shown  presently  that  there  is  a 
difference  in  the  statistics  of  total  number  photo- 
electrons which  are  generated  in  an  interval  T for 
coherent  and  noncoherent  light.  We  define,  therefore, 
the  integrated  sampled  photon  arrival  rate  n^^: 

/t  +T 
m 

N (t)  dt 
e 

^m 


f i (t 


)dt  = 


qM, 


pm 


J 


Sr. 


m 


m 


We  can  relate  n to  the  intensitv  I (t)  by 

m * 

using  the  relationships  presented  above: 


t +T 

CA,  ^ m 


m 


^ r 

hv^  J 


I(t)  dt 


m 


where 


/ 


t +T 
m 


= n J I (t)  dt 

^m 


n = 


hv. 


The  instrumentation  depicted  in  Fig.  1 is  assumed 
to  include  a digital  processor.  This  processor  will 
collect  data  (n^^)  from  M samples  and  compute 
various  statistical  parameters  on  this  sample  space 
as  indicated  below. 


5.  THEORY  OF  OPTICAL  SIGNAL  DETECTION  METHOD 


For  an  understanding  of  the  present  optical 
detection  method  it  is  helpful  to  examine  the  relation- 
ship between  the  statistical  distribution  function  for 
the  photon  arrival  rate  and  the  optical  intensity  I(t). 
For  our  study  we  need  the  probability  p(n,  T)  (i.e., 
ergodic  statistics)  of  emitting  n photo  electrons  in 
time  interval  from  t to  t + T where  T is  fixed. 

It  is  shown  in  Ref.  1 that  for  perfectly  coherent  light 
of  intensity  I(t)  we  have: 


1 

/-t+T  . .1" 

/•t+T  , ,1 

P(n,T,t)  * J5J  n 

/ I(t  )dt  exp  -r 

t J 

) 7^  i(t  )dt  I 

Note  that  this  reduces  to  the  Poisson  distribution  if 
the  mean  number  of  emitted  photo  electrons  is  denoted 
<n> : 


p(n,T)  - ^ (nW)" 


where 


<n>  * nW 


w 


r t+T  , , 

7 lit  )dt 

•'  t 


n = 


h<v> 


? = quantum  efficiency  of  detector 
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<v>  = average  optical  frequency 

h = Planks  constant 

I 

Fcr  noncchaiant  light  I (t  ) is  a stochastic 
process  and  the  photon  counting  distribution  must 
be  averaged  over  all  possible  I(t). 

p(n,  T)  * <(nW')’^  exp  (-nW)> 

n • 


- krf 


(nw)”  exp  (-nW)  p(W)  dw 


We  can  illustrate  the  use  of  this  relationship  by 
finding  p(n,  T)  for  an  ideal  amplitude  stabilized  laser. 
In  such  a case  the  probability  density  function  p(I)  for 
the  intensity  is  given  by: 

P(I)  = 6(1  - <I>) 


where 


<I>  = average  intensity 

and  where  6(-)  is  the  Dirac  delta  function.  The 
integrated  intensity  W is  given  by 

W = IT 


The  photo  detection  equation  becomes: 


= if-  <hW>" 
n! 


-<n  w> 

e 


p(n,  T) 
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<n>*'  _-<n> 

“ “jn-  ® 

where  <n>  = n <IT>  = n <I>T  is  the  average  number  of 
photo  electrons  generating  in  an  interval  of  duration 
T.  Note  that  this  is  the  Poisson  distribution. 

Another  important  special  case  merits  attention. 

We  consider  the  photon  counting  distribution  for  a 
thermal  source.  An  optical  field,  which  is  generated 
by  a thermal  source  composed  of  many  Independent  atomic 
radiators,  consists  of  many  different  frequency  components 
occupying  a fixed  band.  From  the  central  limit  theorem 
it  can  be  concluded  that  the  resulting  optical  field 
corresponds  to  a zero  mean  Gaussian  random  process  in 
scalar  amplitude. 

Assuming  a linearly  polarized  thermally  generated 
optical  field  then  the  complex  scalar  amplitude  can  be 
denoted  U where ; 

U = ^ ^i 

and  where 

= Re  U 

U.  = Im  U 
1 

The  joint  probability  density  function  for  these 
components  are  given  by; 


X J 


2 

where  o is  the  variance  in  the  intensity  I; 

I = ju!^ 

= Var  (jul^)  = ^ 

We  can  write  this  joint  density  function  in 
polar  form  where 


U = |u|e^'^ 

(u|  = /u^  + vl 


U. 

tan  (J)  = 


Letting  (u|  = /T  we  obtain 


P(/I,  <{)) 


(/I)^ 

/l <I> 

n<i>  ® 


This  probability  density  function  is  independent  of 
the  phase  <J)  of  the  scalar  amplitude.  Thus  we  can 
integrate  the  joint  density  function  over  ij)  to  obtain 
the  probability  density  function  for  the  intensity 
alone : 


P(/I) 


2/1 

<I> 


e 


I 

<I> 


li 


I 

I 
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This  function  can  be  further  converted  from  a function 
of  /T  to  a function  of  I by  noting 

dl  = 2 /I  d /I 


Thus  we  obtain 


P(I) 


I 

<I> 


In  our  instrumentation  we  will  be  using  very  fast 
detectors  whose  outputs  will  be  sampled  for  very  short 
intervals  compared  to  the  time  in  which  I(t)  changes 
appreciably.  In  this  case  letting  T be  the  sampling 
interval  we  can  approximate  the  integral  for  W as 
follows 


w 


I(t')  dt 


% IT 


From  this  we  obtain  for  p(n,  T) 


P(n,  T)  = 


/ 


^ 


<n> 

1 + <n> 


1+n 


where  <n> 


n<I>T  is  equal  to  the  average  number  of 


basis  for  our  method  of  detecting  laser  radiation. 

To  enhance  this  notion  it  is  instructive  to  consider 
the  mean  and  variance  for  these  two  distributions.  If 
we  imagine  that  the  light  which  is  incident  upon  our 
detection  apparatus  is  totally  coherent  then  we  can 
denote  the  optical  intensity  I^Ct)  and  the  number  of 
photo  electrons  generated  in  the  sampling  interval  T. 

We  further  presume  that  this  light  is  generated  in  a 
perfect  amplitude  stabilized  laser  such  that  the 
Poisson  arrival  rate  statistics  apply.  Denoting  the 
mean  photon  arrival  rate  the  variance  is  given 

by 


V(n^) 


2 

oo  (n^  - <n^,>)  e n^, 

Z ; <n.> 

n=l  "c*  ^ 


= <n 


> 

c 


This  result  is  a well  known  property  of  the  Poisson 


process . 
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On  the  other  hand,  assuming  that  the  incident  light 

is  totally  noncoherent  we  denote  the  intensity  and 

the  numbers  n^  of  photo  electrons  generated  in  interval  T. 

Denoting  the  average  of  this  number  of  photo-electron 

<n.>  we  obtain  for  the  variance: 

1 


V(n^) 


2 "^i 

(n^  - <n^>)  <n^> 


ni=l 


(1  + <n^>) 


1+n. 


= <n^> (1  + <n^>) 


Note  that  if  we  take  the  ratio  of  variance  to  mean 

for  the  two  cases  we  have  the  basis  for  discrimination 

between  coherent  and  noncoherent  radiation.  We  denote 

this  ratio  and  S.  for  coherent  and  noncoherent 

c 1 

respectively  and  obtain: 


V(n^) 

<n.> 


1 = 0 


V(n.) 

S . = =- 

1 <n.> 

1 


1 + <n^>  - 1 


<n . > 
1 


These  functions  are  sketched  in  Fig.  2.  By  measuring 
the  mean  and  variance  of  the  detected  photo  electrons 
we  can,  in  principle,  distinguish  perfectly  coherent 
from  perfectly  noncoherent  light.  Note  that  it  is 
advantageous  to  pick  T large  enough  such  that 


6.  FORMULATION  OF  SIGNAL  DETECTION  PROBLEM 


In  the  previous  section  of  this  report  it  has  been 
shown  that  the  statistics  of  the  photon  arrival  rate 
are  a function  of  temporal  coherence.  This  influence 


of  coherence  is  expressed  by  the  probability  density 


function  for  the  number  of  photons  n counted  by  a 
photo  detector: 


p(n)  = 


(n 


— m--— 1 ^ ^ 


.m-1 

-n 


<n  >m 
c 


<n^> (<n^>  + m) 


where 


Lm-1 

n 


( ) 


m 


<n . > 
1 


<n  > 
c 


the  Laguerre  polynomial 

number  of  modes  in  the 
optical  field 

average  number  of  noncoherent 
photons  counted  in  an  inter- 
val T 

average  number  of  coherent 
photons  counted  in  T 


The  number  of  modes  m of  the  optical  field  is  deter- 
mined by  the  receiver  optics.  It  is  shown  in  Ref.  1 that 
this  number  of  modes  is  given  by: 


m 


ss"- 

2 2 


T 

T 


X 
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where 

S = detector  area 

= source  area 

R = distance  from  sound  to  detector 

A = optical  wavelength 

T = sampling  interval  in  which 

photon  rate  is 

X = coherence  time  of  source 


In  the  proposed  instrumentation,  the  detector  output  will 

be  sampled  for  a time  T.  The  number  of  photons  in  the 

mth  sampling  interval  of  duration  T is  denoted  n . 

m 

If  this  sampling  is  taken  near  the  nyquist  limit,  then  T 
will  be  of  the  order  of  the  reciprocal  of  the  detector 
bandwidth  6^ : 


T 


• i 


— 6 “9 

A typical  value  for  T is  in  the  range  from  10  to  10 
In  most  applications  the  wavelength  of  the  laser 
being  sought  is  known.  In  this  case  it  is  reasonable  to 
utilize  an  optical  preselector  to  filter  the  light  reaching 
the  detector.  A Fabry  Perot  etalon  is  an  example  of  such 
a preselector. 

Normally  it  is  desirable  to  have  the  preselector  band- 
width as  narrow  as  possible  for  maximum  pre-signal  proces- 
sing signal/noise.  A Fabry  Perot  etalon  is  commercially 
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o 

available  having  a bandwidth  AX  of  about  lA.  The  light 
emerging  from  a filter  of  this  bandwidth  will  have  a 
"coherence  time"  x of  the  order  of 


T 


O 

where  C equals  speed  of  light  and  X wavelength.  A lA 
filter  resonant  at  1 um  will  yield  a coherence  time  of 
about  10  sec. 

It  is  instructive  to  estimate  the  number  of  modes 
in  a filtered  optical  field  for  a typical  set  of  optics 
that  could  be  used  for  detecting  coherent  light  in  a 
noncoherent  background.  The  following  "reasonable"  values 
are  taken  for  the  parameters  in  the  expression  for  m: 


S ~ 
- 


R 

X 


1 cm 
1 m^ 
1 Km 
1 um 


For  this  representative  example  the  number  of  modes  is 
given  by 


m = 100 

In  general  the  number  of  modes  in  an  optical  field  appropriate 
for  the  proposed  signal  detection  scheme  will  be  in  the 

4 

range  100  < m < 10  . However,  at  times,  single  mode  fields 
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! 

i 


will  be  considered  in  the  following  analysis  for  simplicity 

in  interpreting  certain  results. 

The  objective  of  this  investigation  has  been  to  find 

means  of  detecting  coherent  photons  in  the  presence  of 

noncoherent  photons.  This  is  accomplished  by  measuring 

various  moments  of  the  distribution  p(n)  and  by  knowing 

the  influence  of  nonzero  <n^^>  upon  these  moments.  Earlier 

we  showed  that  the  variance  of  p(n)  depends  upon  <n^> 

and  <n  >.  In  doing  so  we  demonstrated  the  variance  for 
c 

the  two  special  cases;  (a)  <n.>  = 0 and  <n  > 5^  0 and 

re 

(b)  <n  > = 0 and  <n.>  ^ 0. 
c 1 

It  is  instructive  to  consider  the  variance  V(n)  for 
the  above  density  function  in  the  case  of  combined  coherent 
and  noncoherent  light  because  this  is  the  applicable  case 
for  the  detection  of  coherent  light  in  a noncoherent  back- 
ground. 

The  variance  of  p(n)  is  most  easily  computed  from 
the  recurrence  relation  for  the  moments  <n  >.  In  Appendix 
A it  is  shown  that 


^ k+1^ 
<n  > 


<n 

c 


+ 


n . ><n^> 
1 


+ 


3<n^> 

3<n^> 


+ 


<n  > 
c 


22V\ 

m / 


3<n^> 

3<n^^> 


where  m is  the  number  of  modes  of  the  optical  field 
being  investigated.  The  coherent  and  noncoherent  photons 


i 

( 


I 

I 


\ 


i-i 


h 
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are  independent  because  they  are  generated  by  independence 
physical  mechanisms.  As  a result  the  first  moment  of  p(n) 
is  given  by 


<n>  = <n  + n.>  = <n  > + <n.> 

Cl  c 1 


From  the  recurrence  relation  the  second  moment  is 
found  to  be: 


<n^> 


- (<n^>  + <n.>)2  . * 2^) 


The  variance  [i.e.,  V(n)]  in  p(n)  is  given  by 


V(n)  * < (n  - <n>) 

2 2 
* <n  > - <n> 


(<n^>  + <n^>)^  + <nj>(l  + 

* ""i")  ^ 


<n  >(1  + 

c V m 


It  is  convenient  to  define  the  ratio  of  average  coherent 
to  noncoherent  photons 


r = 


<n  > 
c 

<n.> 


Using  this  notation  the  variance  V(n)  becomes: 
V(n)  * <n: 


r 

1 + ' 

+ 1 

"l  1 1 

1 + r 

_ m(l+r)_ 

1 + r 

^ m(l+r)J 
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Note  that  this  variance  which  depends  significantly 
upon  r can  serve  as  the  basis  for  detecting  coherent 
light.  Totally  noncoherent  light  corresponds  to  the 
special  case  r = 0.  Totally  coherent  light  corresponds 
to  the  limit  as  r tends  to  infinity.  For  any  given 
signal/noise  (i.e.,  for  any  given  r)  the  detection  of 
coherent  light  involves  distinguishing  between  V(n,r) 
and  V(n,o).  It  is  the  purpose  of  the  following  analysis 
to  consider  the  process. 


6.1  Definition  of  Signal 

In  attempting  to  detect  temporally  coherent  light 
in  a noncoherent  background  from  the  photon  arrival  statis- 
tics, the  received  signal  must  be  sampled.  In  the  proposed 
scheme,  the  nximber  of  photons  is  counted  for  an  interval 
T repeatedly.  Denoting  the  number  of  photons  received  in 

the  mth  sample  M , then  an  estimate  of  the  variance  V 

m 

of  n can  be  obtained 


V(n) 


N - 1 


N 

Z 

m=l 


m 


<n> ) 


2 


where 


<n> 


n 


m 


i.e. , <n>  is  an  estimate  of  the  mean  number  of  photons  in 
an  interval  T and  where  N is  the  number  of  samples 


f: 

li-i 


! 


taken. 
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For  mathematical  convenience  we  assume  that  N is 
sufficiently  large  that  the  estimate  of  <n>  is  essentially 
an  accurate  measure  of  <n>: 


<n> 


<n> 


Lim 

N-voo 


<n> 


It  is  inadequate  to  consider  the  variance  of  n for  our 
purposes.  However  the  relationship  between  variance  and 
mean  of  n can  serve  as  the  basis  for  detection  of  coherent 
light. 

In  Section  5 of  this  report  we  introduced  a definition 
of  signal: 


S » ^ 

<n> 


1 


We  use  a similar  definition  here  for  S,  which  can  be 
written  in  the  following  form: 


S 


V(n) 

<n> 


<n>  1 + 2r 
(1+r)^ 

It  is  important  to  recognize  that  the  proposed 
instrumentation  scheme  for  detecting  coherent  light  is 
capable  of  only  estimating  S.  That  is  denoting  this 
estimate  S,  the  proposed  instrumentation  can  perform 

A 

computations  leading  to  a value  for  S 


xl 
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= _ 


The  quantity  denoted  S,  which  we  have  called  signal,  is 
the  expected  value  of  this  estimate: 


S = <S>  = 


Lim  S 

N-vco 


6.2  Post  Signal  Processing  Noise 

For  any  finite  sample  space  (i.e.,  finite  N)  the 
/\ 

estimate  of  S will  be  in  error  by  an  amount  denoted  e 


e = S - S 


where  e is  a random  variable.  A convenient  measure  of 

A 

the  size  of  the  error  e is  the  variance  in  S,  i.e., 


V(S)  = < (S  - S)^> 


= <e  > 


Tv(n)  , V(n)  . ,T 
fvln)  V(n)]^. 

[ZT  - — J 


In  principle  we  could  expand  the  above  expression  in 

A A 

terms  of  the  variance  of  V and  <n>.  However,  the  compu- 
tations are  unwieldy  and  will  yield  results  that  are  not 
readily  interpreted.  Fortunately  for  relatively  large 
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sample  space  (i.e.,  large  N)  the  variance  <n>  is  small 
compared  with  the  variance  in  V: 

v[<n>]  = <<  v[v(n)] 

The  validity  of  the  above  inequality  can  be  established 
by  numerical  exaimple.  In  this  case  the  following  approxi- 
mation can  be  made: 

<n>  s:  <n> 

Substituting  this  approximation  into  the  expression  for 

A 

V(S)  yields 


V(S)  “ — v[v{n)] 

<n> 


It  is  shown  in  Ref.  3 that  v[v(n)3  is  given  approxi- 
mately by 


N 


2u 


2 

2 


where  and  ^2  the  fourth  and  second  central 

moments  of  p(n)  respectively: 

= < (n  - <n>) 

“ 4 

= E (n  - <n>)  p(n) 

n=l 

2 

U2  = < (n  - <n>)  > 
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U,  = E (n  - <n>)^  p(n) 

^ n=l 

It  is  possible  to  evaluate  and  U2  from  the 

moments  <n  > (k  = 1,  2,  3,  4)  using  the  recurrence  rela- 
tionship that  was  developed  in  Appendix  A.  In  Appendix  B 
it  is  shown  that  the  variance  in  the  signal  is  given  by: 


V 


N<n; 


* * .H 

[i^(' 


2<n>  \ 

<n>  ( 

1 + ^ 

m(l+r)y 

1 + r V 

m(l+r)y 

2<n>  \ 

, 3<n>  1 

m(l+r)7 

m(l+r) 

m(l+r) 

r) 


In  characterizing  measurement  errors  it  is  common  practice  to 
specify  the  standard  deviation  a (i.e.,  square  root  of  variance): 


a(S) 


/ 


V(S) 


If  the  proposed  signal  detection  method  were  implemented  with  ideal 
instrumentation  (i.e.,  zero  instrument  noise),  then  the  noise  asso- 
ciated with  this  scheme  results  from  the  finite  sample  space  from 
which  S is  estimated.  A convenient  measure  of  this  noise  is  a(S) 
defined  above. 

A 

■N  = a(S) 
o 


It  is  shown  above  that  S is  given  by 


<n>  1 + 2r 
(1+r)^ 


Figure  3 is  a sketch  of  signal  S vs  <n>.  Note  that 

S(<n>)  is  a linear  function  of  <n>  having  a slope 
1 + 2r 

=•  for  any  given  ratio  r and  number  of  modes  m. 

M(l+r) 

The  essential  features  of  our  signal  detection 
problem  can  be  understood  with  respect  to  Fig.  3.  The 
instrumentation  must  distinguish  between  a signal  having 
r identically  zero  (i.e.,  no  coherent  light)  and  light 
having  small  but  nonzero  r.  The  slope  of  the  S(<n>) 
curves  is  a decreasing  function  of  r having  the  maximum 
value  at  r = 0 (i.e.,  zero  signal)  and  having  zero  slope 

for  r-*-«  (i.e.,  totally  coherent  optical  field). 


Figure  3 
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The  signal  detection  problem  involves  a measurement 
of  r in  order  to  determine  nonzero  r.  If  it  weren't 
for  noise,  the  signal  detection  problem  would  be  straight- 
forward. A suitable  signal  detection  instrument  would 
measure  <n>  and  then  S<n>  for  the  given  (and  fixed) 
m could  be  computed.  From  this  r could  be  determined 
by  inverting  the  expression  for  S(<n>,  r) . 

Unfortunately  there  is  considerable  noise  associated 
with  the  proposed  detection  scheme,  which  places  a lower 
limit  upon  detectable  r.  As  mentioned  above,  there  is 
optical  receiver  noise,  which,  for  convenience,  is  being 
neglected.  The  fundamentally  unavoidable  noise  associated 
with  the  proposed  signal  detection  method  results  from 
the  finite  sample  estimate  of  S as  explained  above. 

This  noise  explains  the  uncertainty  in  the  measurement  of 
S thereby  placing  a limit  on  the  minimum  r which  can 
be  measured. 

6.3  Excimple  Signal  Detection  Method 

For  an  understanding  of  the  limits  placed  upon  the 
minimum  size  of  r which  can  be  measured  it  is  helpful 
to  refer  to  Fig.  4.  This  Fig.  4 is  a sketch  of  the 
probability  density  function  for  S.  This  function  tends 
to  express  the  distinction  of  measurement  results  for 
S(r)  for  any  given  true  S(r).  Here  we  have  depicted 

/A 

the  conditional  probability  for  S under  the  two 


I 


S(0) 

Figure  4 

hypotheses  r = 0 and  r > 0.  The  variance  for  these  two 
conditional  density  functions  are  v[s(0)]  and  v[s(r)] 
respectively.  Recall  that  we  have  shown  these  variances 
to  be  inversely  proportional  to  the  size  of  the  sample 
space  N. 

6.4  Feasibility  of  Proposed  Methods 

In  the  following  subsection  of  this  report  the 
feasibility  of  this  proposed  signal  detection  method  is 
considered  with  respect  to  certain  practical  limitations. 
In  considering  this  feasibility,  it  is  helpful  to  consider 
a specific  example  signal  detection  method. 

In  one  very  traditional  implementation  of  a signal 


A 


1| 


detection  scheme  a simple  threshold  decision  is  reached. 

The  decision  will  be  reached  that  signal  is  present  whenever 


< CO 
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where  is  the  threshold  signal-to-noise  ratio  for 

reaching  the  decision  that  the  signal  is  present.  Referring 
to  Fig.  4 this  criterion  for  signal  present  is  equivalent 
to  the  condition 


S(r)  < 

where  = S(r^).  These  conditions  are  equivalent 

because  S is  a monotonic  nonincreasing  function  of  r. 

/s.  /N 

In  the  proposed  implementation  both  <n>  and  S 
would  be  obtained  from  N samples  2,...,  N)  . 

In  our  above  analysis  it  has  been  assumed  that  N is 

/s 

sufficiently  large  that  <n>  « <n>.  With  known  <n> 
and  m (i.e.,  m fixed  by  configuration),  it  is  possible 
to  compute  S(0) 


S(0) 


<n> 

m 


The  minimum  desirable  detectable  signal  level  r . 

^ min 

is  selected,  and  S . = S (r  . ) is  then  determined. 

min  min 

The  threshold  level  will  be  somewhere  between  these 

two  values  of  S. 


S . 
min 


< < S(0) 


depending  upon  the  choice  of  acceptable  errors. 

As  in  all  such  signal  detection  methods  there  are 
two  classes  of  error  which  can  be  made:  (1)  failure  to 
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detect  the  signal  when  it  is  present  and,  (2)  deciding 

the  signal  is  present  when  there  is  none.  The  probability 

of  an  error  of  the  first  (i-e.,  P ,)  is  given  bv 

ei 

'■el  ' /" 

^t 

and  the  probability  of  the  second  type  is  (i.e.,  ' 

/S  /V  /V 

„ f2(S)dS 

The  relative  size  of  these  two  error  probabilities  depends 
upon  the  function  Pj^(S)  and  ^2^^^  upon  the  choice 

of  S^.  Selection  of  an  optimum  value  for  is  well 

covered  in  the  general  theory  of  signal  detectability  and 
is  not  considered  here. 

Although  this  investigation  is  not  a study  of  signal 
detection  theory,  we  can  consider  a simplified  and  interes- 
ting special  case  for  illustrative  purposes.  We  assume  that 
Pj^  and  P2  have  identical  shapes  (i.e.,  identical  central 
moments)  but  different  mean  values  [i.e.,  S(r)  5^  S(0)]. 

We  further  assume  for  mathematical  convenience  that  it  is 
desirable  for  p^^^  = ^e2*  this  is  the  criterion  for 

errors,  then  should  be  halfway  between  and 

S(0)  : 


S 


t 


S(0)  + s 


min 


2 
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Furthermore  the  total  probability  of  error  in  our 

simplified  example  is  the  following  sum: 


= P 


el 


e2 


/o  ^ ^ 

^ p,(S)dS  + 
0 ^ 


//\  /V 

p,  (S)dS 


For  any  given  pair  of  probability  density  functions 

p^  and  P2  (e.g.,  Gaussian)  there  is  a unique  relationship 

between  P_/  S.  and  the  moments  of  these  functions.  For 
E t 

example  P„  is  a function  of  S.  and  a.  o_  (i.e.): 

Pg  = f[S^,  a^iS),  a2(S)] 


where  and  02  are  the  standard  deviations  of  pj^ 

and  P2  respectively. 

In  a typical  signal  detection  problem  a positive 

real  value  between  zero  and  unity  is  assigned  to  P„ 

£« 

representing  an  acceptable  level  of  error.  Having  previously 

determined  the  minimum  acceptable  ratio  of  signal  to  noise 

(r  . ) and  S^,  we  can  determine  the  maximum  values  for 
min  t 

o.  and  a_  from  the  allowed  P„.  From  these  maximum 

1.  Z £« 

standard  deviations  we  can  determine  the  required  size  of  the 

scimple  space  (N)  . The  feasibility  of  the  proposed  signal 

detection  scheme  can  be  assessed  from  this  value  for  N. 
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There  are  many  other  approaches  for  assessing  feasi- 
bility based  upon  other  sets  of  constraints.  For  example, 
the  sample  space  size  can  be  fixed  by  the  sampling  rate  of 
the  receiver,  by  the  duration  of  the  transmitter  coherent 
light  signal  and  by  the  computational  capacity  of  the 
signal  processing  unit.  However,  we  have  chosen  the  above 
criteria  for  feasibility  for  illustrative  purposes. 

Understanding  of  this  feasibility  analysis  is, 
perhaps,  enhanced  by  specific  numerical  examples.  In  these 
examples  we  select  a fixed  ^ specific  optical  power 

level  p^  wavelength  X bandwidth  Av  and  detector  band- 
width 8^.  From  these  we  compute  the  required  sample  size 

N frOTi  which  S is  estimated.  We  select  r„.  = 0.1 

min 

which  is  equivalent  to  -10  dB  (power  ratio)  signal/noise. 

In  the  following  examples  the  probability  density  functions 

Pj^  and  P2  will  be  approximated  by  Gaussian  functions  for 

S ^ 0.  This  approximation  is  reasonable  provided  the 

variances  of  p^  and  p^  are  sufficiently  small  or 

equivalently  provided  N is  sufficiently  large.  In  our 

first  example  the  total  received  optical  power  is  taken  to 

be  10~^^  watts  at  an  average  photon  energy  <hv>  of 

10  Joule  (i.e.,  near  infrared).  We  assume  that  the 

“8 

sampling  interval  is  T = 10  sec  which  is  equivalent  to 
a nyquist  seunpling  rate  of  a detector  having  a bandwidth 
of  the  order  of  50  MHz.  The  average  number  of  photons 
received  in  time  T is  given  by 


! 
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<n> 


Pt 

<nv> 


= 100 

In  attempting  to  assess  feasibility  of  this  method  with 
the  following  special  cases  we  will  refer  to  "large"  and 
"small"  value  for  <n> . By  such  terms  we  imply  large  or 
small  compared  with  <n>  = 100  photons/sampling  interval. 
The  implication  of  the  magnitude  for  <n>  involved  in  any 
particular  case  with  regard  to  optical  power  level  and 
bandwidth  can  be  determined  from  a numerical  estimate 
similar  to  the  preceding  example. 


6.5  Special  Case  - Single  Mode  Fields 

Although  it  will  presently  be  shown  that  it  is  not 
a physically  realistic  example,  single  mode  (m  * 1)  fields 
are  most  easily  interpreted.  The  variance  V(S)  must  be 
determined  for  our  analysis.  For  signle  mode  fields  the 
leading  term  in  V(S)  for  large  <n>  is 


2<n>  ^ 

1 + r 


The  standard  deviation  and  02  are  given  by 


a 


1 


<n> 

1 + r 


r'W 


/ v[S(r)] 


Note  that  for  r - 0.1  these  two  standard  deviations  are 
approximately  equal,  i.e.. 


01  » 02 


The  criterion  for  determining  N can  be  understood 
with  reference  to  Fig.  5 which  is  a sketch  of  the  two 
probability  density  functions  p,  and  p,. 


Figure  5 


Recall  that  the  total  probability  of  error  P_  is  given 

t* 


/S » A ^ ^ ^ ^ 

^ p^(S)dS  + J p,  (S] 
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Pj^(S)dS 

For  the  assumed  Gaussian  function  P_  is  given  approxi- 

Cl 

mately  by 


45 


Solving  the  above  equation  for  the  maximum  N t^niin^ 
yields 


N 


(4.8)^ 

4 


^rain 


25  X lo"* 


The  feasibility  of  the  proposed  method  is  influenced 

by  the  memory  capacity  of  the  signal  processing  and  the 

duration  of  the  transmitted  signal  from  which  the  N scunples 

“8 

are  taken.  If  the  s«imples  are  taken  in  intervals  T * 10 
sec,  then  the  minimxim  signal  duration  for  x for  all  N 
samples  must  be 


T » lO"®  N 
* 2.5  msec 

In  most  practical  cases  a signal  duration  of  2.5  msec 
will  not  be  limiting.  However,  a required  memory  capacity 
to  store  250,000  data  points  is  totally  unrealistic.  Thus 
it  appears  that  the  implementation  of  the  proposed  scheme 
using  a computer  or  special  purpose  processor  is  not 
feasible. 

On  the  other  hand,  a minimum  signal/noise  of  0.4 
requires  a memory  capacity  of  only  about  N ~ 1000  samples 
which  is  clearly  within  the  capability  of  a small  digital 
processor. 

Although  our  example  feasibility  study  is  not  physically 
realistic,  it  does  demonstrate  a few  of  the  fundamental 


V 


I 
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compromises  between  system  complexity  and  performance  of 
the  proposed  method.  The  sensitivity  of  the  size  of  the 
sample  space  upon  the  required  memory  capacity  is  clearly 
demonstrated.  The  following  numerical  example  further 
expands  the  understanding  of  the  influence  of  system  para- 
meters upon  performance. 


6.6  Second  Special  Case  Feasibility  Stud 


Our  second  special  case  estimate  of  the  feasibility 
of  the  proposed  signal  detection  method  considers  multimode 
fields.  In  this  second  special  case  we  once  again  assume 
relatively  large  number  of  received  photons  (<n>  = 100)  in 
the  sampling  interval  T. 

The  number  of  modes  m associated  with  the  optical 
field  can  be  computed  from  the  configuration  depicted  in 
Fig.  6.  It  has  been  shown  above  that  m is  given  by 


SS^  Av 

“2  5 
R 


We  assume  the  following  reasonable  values  for  the  parameters 
of  the  expressions  for  m: 


= 0 . 1 m'" 

1 2 
= 10  m 


= 1 Km 


= 1 ym 


Source 

Area 


Figure  6 


Substituting  these  numerical  values  into  the  equation 
for  n yields: 


m = 


10”^  X lo'^  X 10^ 
10®  X 10-^^ 


For  the  above  values  of  <n>  and  m the  leading 
term  for  V(S)  is  given  by 
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The  standard  deviation  for  and  P2  is  given  by 


o = 


/N 


For  a total  = 0.01  we  have  once  again  (assuming 

Gaussian  functions) : 


S(0)  - S(r)  = 4.8  o 


<n>  ( r y _5_ 

V*  ■ /if 


Solving  for  N yields; 


Taking  r„.  = 3 and  using  the  above  determined  values 

min 

for  m and  <n>  yields: 


N = 7500 

This  is  a reasonable  sample  size  for  a small  general 

purpose  computer  or  processor.  If  the  memory  capacity  is 

increased  to  16  K,  then  r . is  reduced  to  only 

min 


min  r z 2.3 


Thus  it  appears  that  the  proposed  method  is  applicable 
for  signal/noise  of  about  6 dB. 
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6.7  Third  Special  Case  - Small  <n> 

The  third  special  case  for  which  the  feasibility 
of  the  proposed  method  is  investigated  is  for  small  <n> 
(small  compared  with  unity) . The  leading  term  in  the 
expression  for  V(S)  in  this  case  is: 


V(S) 


<n>-»-0 


1 

<n>  H 


{<n>} 


1 

<n>N 

The  standard  deviation  for  this  case  is  given  by: 

a = 

/<n>N 


The  criterion  of  a 0.01  total  Pj^  is  then; 


s<0)  - 


= 2 a 


/<n>N 


Solving  for  N yields; 


N = 


4 m" 
<n>‘ 


In  interpreting  this  equation  it  is  sensible  to  fix 

a reasonable  value  for  N„  and  to  find  the  value  for 

max 


lasgsiBKKnaoBJCB 
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<n>.  We  have  shovm  that  for  a typical  configuration  m 2 10 
A reasonable  upper  bound  for  N is  of  the  order  of  16  x 10 
If  we  arbitrarily  select  r = ^ , we  find 


_ 1 

<n>  = (20)  ^ = 0.3 


This  corresponds  to  an  optical  power  of 


P = 


<n> 


hv 


-20 

10  X 0.3 


10 


-8 


-12 

0.3  X 10  watts 


6.8  Fourth  Special  Case  - Optimiim  <n> 

The  next  special  case  which  we  consider  involves  the 
selection  of  <n>  for  minimxim  variance  in  our  estimate  of 

/S 

signal  S.  For  a given  optical  power  level  we  can  choose 
<n>  by  selecting  the  seimpling  interval  T (neglecting 
physical  constraints) . 


<n> 


P 

<hv> 


T 


It  can  be  shown  from  an  upper  bound  for  V(S)  that 
there  is  an  optimum  <n>  which  minimizes  the  variance  and 
hence  the  standard  deviation  a. 

As  we  have  already  shown  the  standard  deviation  repre- 
sents the  post  signal  processing  noise.  In  the  following 
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analysis  we  compute  the  average  number  of  photons  which 
are  counted  in  the  sampling  interval  T for  which  V(S) 
is  minimized.  This  simultaneously  minimizes  the  standard 
deviation  and  noise.  Thus  the  present  special  case  can 
be  considered  the  minimum  noise  case. 

In  Appendix  C the  following  approximation  is  estab- 
lished  for  V(S)  for  single  mode  fields; 

V(S)  < 2 

Note  that  this  function  asymptotically  approaches  infinity 
for  <n>  -*•  0 and  for  <n>  ^ Hence  there  must  be  a 

/N 

positive  real  <n>  which  gives  minimum  V(S)  and  hence 
minimum  noise. 

A 

Letting  A = N V(S)  and  a = <n>  we  can  find  the 
value  of  a which  minimizes  A by  differentiating; 

It  - 0 - 4 . 2a)  (e  . i)  - 

= 24  a^  + 2a  - 1 

= (6a  - 1)  (4a  + 1) 

The  minimum  for  A occurs  for 

1 

^ = 6 

which  yields 


min  V (S) 


N 
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4 

3 


(12) 


N 


ii 

3N 


It  is  further  shovm  in  Appendix  C that  for  multimode 

/s 

fields  the  variance  V(S)  is  given  approximately  by 


V(S) 


(l  + 

1 4-  ^ 4-  C 

r + -L.‘ 

(1  ^ m j 

X T — T o 

m 

^ ^ <n> 

It  is  also  shown  there  that  this  variance  is  minimized  for 


a = 


+ 16  bm 


m 


where 


a = value  of  <n>  which  minimizes  V(S) 


b = 


1 + ^ + 6r 
m 


Earlier  we  have  shown  that  a practical  realization  of  the 
proposed  signal  detection  method  involves  m of  the  order 
of  100  to  1000.  For  m greater  than  about  100  a(m) 
is  given  approximately  by 


a (m) 


2 11  + 


6r 


--^11 


Figure  7 is  a graph  of  a(m)  for  a few  representa- 


tive values 


100 


10 


I (m) 


1 


Notice  that  as  the  signal/noise  increases  for  any  given  m, 
the  required  average  photon  count  (<n>)  decreases  for  the 
minimum  noise  case.  Notice  further  that  in  the  normal 
operating  region  (100  < m < 1000  and  r < 1)  the  average 
number  of  photons  in  a sampling  interval  T is  about  10. 

/s 

For  m > 100  (the  typical  case)  the  variance  V(S) 


is  given  by 
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V(S) 


1 6r 
N 


The  total  probability  of  error  P_  will  be  less  than  or 

L 

equal  to  0.01  provided 


The  required  sample  space  size  for  this  case  is  given 
by 


N 


25(H-6r) 

IJIKZ 

2 ll+6r 


3 

= 50(1  + 6r)^ 


If  the  memory  capacity  is  limited  to  16  K,  then  the 
minimum  r is  given  by  about 


min  r ::  1.4 


This  shows  an  improvement  of  about  a factor  of  2 over  the 
minimum  detectable  signal/noise  in  the  large  <n>  special 
case.  Thus  there  are  benefits  to  operating  at  optimum  <n>. 

It  is  instructive  to  consider  that  the  mechanism  by 
which  <n>  is  selected  for  any  given  received  optical  power 
P is  made  by  selecting  T.  The  relationship  between  <n> , 
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P and  T is  given  by 


<n> 


PT 

h'-'v> 


where  <v>  is  the  average  optical  frequency.  This  fre- 
quency is  given  by  the  center  frequency  of  the  received 
optical  band.  Typically  this  will  be  the  band  center  of 
any  optical  preselector  filter  which  precedes  the  detector. 

For  minimum  noise  and  large  m fields  the  average 
number  of  photons  received  in  interval  T is  given  approxi- 
mately by 

” h<v>  “ 2 Vl  + 6r 

We  are  most  interested  in  detecting  temporally  coherent 
light  at  low  signal  levels  and  at  low  signal/noise  (small 
r) . Thus  for  any  given  power  level  P,  the  optimum  choice 
of  T is 


T 


h<v> 

2P 


/m 


Of  course,  this  result  was  obtained  without  regard 
to  any  constraints  imposed  by  the  detector  bandwidth  or 
the  speed  of  the  sampling  electronics.  In  any  practical 
implementation,  operation  at  the  minimum  noise  condition 
may  not  be  feasible. 


7.  SUMMARY  OF  RESULTS 


This  report  has  described  studies  of  a method  of 
detecting  temporally  coherent  light  in  a noncoherent  back- 
ground. It  has  been  shown  that  the  arrival  rate  statistics 
of  photons  are  different  for  coherent  and  noncoherent 
light.  In  this  study  the  probability  density  function 
p(n)  for  the  member  n of  photons  which  arrive  in  an 
interval  T is  a function  o^.  the  ratio  r of  the  average 
coherent  to  noncoherent  photons. 

In  addition  it  has  been  shown  that  the  various 
moments  of  p(n)  depend  upon  r.  From  this  observation 
a post  signal  processing  signal  S has  been  defined: 


S 


V(n)  _ , 
<n> 


where 


V(n)  = <(n  - <n>)^> 

and  where  < > indicates  expected  value.  It  has  been 
shown  that  S is  a function  of  r.  Moreover  it  is  shown 
that  in  principle  it  is  possible  to  distinguish  the 
presence  of  temporally  coherent  light  in  a noncoherent 
background  (r  > 0)  from  purely  noncoherent  light  (r  = 0) . 

A perfect  measurement  of  r would  require  a zero 
error  measurement  of  S.  Unfortunately  it  is  impossible 
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to  measure  S without  error.  Rather,  S can  only  be 
estimated  from  a finite  s^unple  space  of  measurements  of 
the  number  of  photons.  This  estimate  of  S is  denoted 
S: 


The  post  signal  processing  noise  associated 

with  a measurement  of  S is  given  by 


N 


o 


o(S) 


The  required  size  of  the  sample  space  N for  detection 
of  coherent  light  at  a given  r is  determined  for  a given 
allowable  error  rate.  The  feasibility  of  this  proposed 
signal  conditioning  is  assessed  by  determining  the  physical 
parameters  involved  in  the  required  sample  space  size.  The 
limitations  on  the  feasible  Scimple  space  size  include;  the 
memory  capacity  and  computation  speed  of  the  processor; 
the  sampling  rate,  detector  bandwidth,  and  coherent  signal 
duration. 

It  has  been  shown  that  for  small  r unimodal  fields 
the  required  sample  space  size  is  given  by 


If  the  memory  capacity  is  to  be  limited  to  about  16  K, 
then  the  minimum  r which  can  be  detected  is  given  by 

min  r ~ 0.11 


Next  this  study  considered  multimode  fields  in  which 
m > 100.  It  was  shown  that  relatively  large  <n>  (about 
100)  the  required  sample  space  size  is  given  by 
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N 


As  an  illustrative  example,  it  was  shown  that  for  min 
r ~ 3 the  required  sample  space  si2e  is: 


N = 7500 


Although  a memory  capacity  of  7500  is  relatively 
small,  increasing  N to  16  K only  reduces  min  r to 
about 


min  r I s:  2 . 3 

I N = 16  K 

Thus  it  appears  that  the  present  method  when  applied  to 
multimode  fields  having  large  <n>  is  limited  to  a signal/ 
noise  of  about  7 dB. 

On  the  other  hand  it  was  shown  that  there  is  an 
optimum  <n>  which  minimizes  noise.  For  m > 1000  that 
value  of  <n>  is  given  by 

1 r~m 

” ~ 2‘yi  + 6r 

For  this  case  it  was  shown  that  the  required  sample  space  is 
given  by 


N = 50 


If  the  memory  capacity  is  limited  to  about  16  K,  then  the 
minimum  detectable  r (for  P^,  = 0.01)  is  given  by 
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min  r ::  1 . 4 

This  is  equivalent  to  about  1 dB  signal/noise  which  is  ’ » 

i' 

roughly  a 6 dB  improvement  over  the  large  <n>  case. 

Thus  it  has  been  shown  that  there  are  benefits  to  opera- 
tion at  optimum  <n>.  ^ 


I 

I 

L 

r< 

n 
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APPENDIX  A 


It  is  the  purpose  of  this  appendix  to  derive  the 
recurrence  relations  for  the  moments  of  the  probability 


density  function  p(n): 


p(n)  = 


1 

(n+m-l)  ! yy  ^^m/ 


l;;'  ^(Z) 


(A.l) 


where 


<n  > 
c 

<n.> 


y (y+m) 

* number  of  modes 


= Laguerre  polynomial 


For  the  purpose  of  establishing  recurrence  relations 


it  is  helpful  to  compute  the  following  derivatives: 


(*> 


feL-'x 


L“  ‘(Z) 


(A. 2) 


mx  1 1+- 


y(l+^)  (n+m-1) 


/ X \ 

v^mAm-1 ' 


(A.  3) 

..\2  m 


L‘"  " (Z) 
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I 

n=l 


= X 


Z [n  - (x+y)]p(n) 
n=l 


k+ 1 k 

<n  > - (x+y)<n  > 


This  can  be  rewritten  in  a more  easily  recognized  form: 


<n  > (x+y) 


Equation  A. 9 is  the  recurrence  relation  which  we  are 
seeking. 


(A. 8) 

(< 
• " 

] 

i 


0 


(A. 9) 


APPENDIX  B 


It  is  the  purpose  of  this  appendix  to  derive  an  expres- 

/s 

sion  for  V(S).  This  represents  the  mean  squared  noise 
associated  with  the  proposed  signal  processing  method.  In 
the  text  it  is  shown  that 


V(S) 


N<n>^ 


(B.l) 


where 

4 

= < (n  - <n>)  > 

2 

^2  = < (n  - <n>)  > 

These  last  central  moments  of  p(n)  can  be  computed  from 
the  first  four  moments  of  p(n)  using  the  recurrence  rela- 
tion of  Appendix  A; 

4 3 2 2 3 4 

* <n > - 4n  <n>  + 6n  <n>  - 4n<n>  + <n>  > 

4 3 2 2 4 

a <n  > - 4<n><n  > + 6<n>  <n  > - 3<n> 


and 


y 


2 


<n  > - <n> 


2 


In  Appendix  A it  is  shown  that 


(B.2) 
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where 


V 

(x+y)<n  > 


3<n^> 

3x 


+ 


3<n^> 

3y 


<Hc> 


<n.  > 
1 


<nc  + n.> 


m = 


X 

y 

X + y = <n> 

number  of  modes 


The  required  moments  can  be  found  in  ascending  order  from 

2 

Eq.  B.3.  First  we  compute  <n  >: 


n^  = (x+y)2  + 

For  the  purpose  of  computing  <n^>  we  find  the  derivatives: 

^ = 2(x-^y)  ^ 

^ = 2 (xH-y)  H-  2|  * 

We  then  obtain 


(B.3) 


(B.4) 


(B.5) 


(B.6) 


(B.7) 
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The  derivatives  of  Eq.  B.7  are  needed  to  find  <n  >; 


From  B.l,  B.8,  and  B.9  we  obtain 


68 


We  obtain  by  substituting  Eq.  B.IO,  B.7  and  B.4 

into  Eq.  B.2: 


Similarly  we  find  ^25 

“2  = 5' (1*5)  * 


(B.ll) 


(B.12) 

X 


The  post  signal  processing  noise  involves  the  binomial 
2 

- 2u2  which  is  given  by 

- 2ul  = 

+ y (1*^)  + X 

The  noise  associated  with  the  proposed  signal  detection 
method  is  most  easily  interpreted  if  we  eliminate  x and 
y from  B.13  in  favor  of  system  variables.  We  define  the 
input  signal  to  noise  ratio  r as 


I m ) m I m I 


I <|-n  -Hi  III-'  II 


ii  nrfi  ^ ^ . 
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t 


r = 


X 

y 


<n  > 
c 

<Hi> 


r 

Noting  that  <n>  = x + y we  find 


r<n> 

^ = r+T 


1 + r 


Substituting  these  expressions  for  x and  y into  Eq. 
B.13  we  can  obtain  v(S)  from  Eq.  B.l: 


APPENDIX  C 


It  is  the  purpose  of  this  appendix  to  show  that  there 

is  an  optimum  average  number  of  received  photon  <n>  in 

the  sampling  interval  which  minimize  the  post  signal 

orocessing  noise  N . It  is  shown  in  the  text  that 
^ o 

N^  = ^v7s) 


and  as  V(S)  is  monotomic  nondecreasing  function  of 

A. 

V(S)  ^ 0,  then  the  value  of  <n>  v;hich  minimizes  V(S) 

also  minimizes  N . 

o 

/N 

The  expression  derived  for  V(S)  in  Appendix  B is 
complicated  and  somewhat  difficult  to  interpret.  On  the 
other  hand,  this  expression  has  an  upper  bound  which  very 
closely  approximates  V(S)  for  a large  number  of  modes 

m.  Recall  that  any  practical  implementation  of  this  method 

. 2 3 

involves  n of  the  order  of  10  to  10  . Thus  it  is 

meaningful  to  consider  the  upoer  bound  for  the  purpose 

of  finding  optimum  <n>. 

In  Appendix  B it  is  shown  that  V(S)  is  given  by 
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where  we  have  used  the  inequality 


, <n>  , 2<n> 

1 + < 1 + - - 
m m 


For  convenience  we  define  the  uoper  bound  for 
A/N: 


V(S)  < g 


where 

^ " 2(1  + (1  + ~ + 6r  + 

Letting  the  value  of  <n>  which  ninir'.ices  A 
a we  obtain 


V(S)  as 


be  denoted 
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(C.6) 


The  value  of  which  minimizes  A can  be  found  bv  differen- 


tiating; 


s • • 


■ »(‘  • *)  ( 


1 + 6r  + ^ + ± 
m a 


(C.7) 


1 + 2 


where 


(C.8) 


1 + ^ + 6r 
in 


(C.9) 


Solving  for  a yields 


-1  + /I  + 4bo 


(C.IO) 


For  large  m (the  most  important  practical  case)  a (m) 
anoroaches 


a (m) 


1 J ^ 

2 f 1 + 6r 


The  result  is  useful  for  interpreting  the  optimum  value 
of  <n>. 


